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Abstract. We study resonant tunneling through a layered medium with a negative 
index medium (NIM) slab as a constituent layer. We demonstrate large delays in 
transmission mediated by the surface and the guided modes of the structure with low 
losses. We show how important it is to include NIM dispersion for correct assessment 
of the nature and magnitude of the delay. We also point out the role of NIM absorption 
for the feasibility of such compact delay devices. 
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1. Introduction 

In recent years there has been a great deal of interest in negative index materials (NIMs), 
which can exhibit exotic properties [H [21 El H]. The rich physics of these materials 
(not occuring in nature) were discussed theoretically by Veselago [TJ, though their 
potential application for beating the diffraction limit [5] and experimental realization 
in the microwave range [6] triggered off the vast current activities. Now such 
materials (also known as metamaterials or left handed materials) have been fabricated 
in other wavelength domains right upto the visible range [31 0]. Applications of 
these metamaterials now range from super-lensing and super-resolution p2, [7] to lasing 
spasers [8] and optical nanocircuits |9j ; invisibility cloaks to electromagnetically induced 
transparency and slow or stopped light [TOl ITT] , etc. The real challenge now is to 
fabricate a low loss metamaterial with high figure of merit (FOM=-^e(n)//m(n) where 
n is the refractive index of the metamaterial) and to broaden the frequency domain where 
both the permittivity and permeability are negative. It is also important to push the 
domain to even higher frequencies maintaining low losses. 

In the context of slow light using NIM, the role of dispersion cannot be 
underestimated. As pointed out by Veselago, a lossless NIM is essentially dispersive 
[!J. The important role played by dispersion has been felt by others in the context of 
waveguides [T2] and especially with reference to cavity QED applications [331 HH fT5] . 
In the context of any finite structure, the frequency dependence of its important 
characteristics (e.g., transmission through it) emerges from two sources [IE]. The first 
is due to the material dispersion of its constituents, while the second is due to the 
boundary conditions. For example, the transmission through an empty Fabry-Perot 
(FP) cavity, or the cavity filled with a dispersive material, can be quite different. These 
characteristics can be quite different as compared to the case of propagation through 
the bulk dispersive sample. We refer to the first as material or intrinsic dispersion, while 
the second is labelled as structural dispersion. In this paper our goal is to bring out 
the salient features of how these sources of dispersion affect the transmission through 
a NIM guide in a resonant tunnelling (RT) geometry [TT] [18] . Unfortunately, many of 
the papers on the NIM waveguides do not address the issue of material dispersion for 
calculating the overall dispersion features [19J. In many others the system is investigated 
at a particular wavelength, concentrating on structural dispersion only, avoiding thus 
all material dispersion related issues [20] . 

Our choice of the RT geometry is motivated by several facts. It was shown recently 
that light can be slowed down using a gap plasmon guide in RT configuration [IS] . Very 
recently a multilayered metal/dielectric structure in the same configuration was studied 
experimentally to show enhanced transmission mediated by the modes of the structure 
[17] . We show that a NIM guide in RT configuration can lead to significant delays 
provided the losses are low. Our calculations are carried out with experimental data for 
the permittivity and the permeability of the NIM [21], albeit with the approximation 
of low losses. We also assumed a homogeneous and isotropic character for the NIM. We 
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studied both sin /cos and the sinh /cosh types of modes of the NIM guide [TH]. We refer 
to the former as guided modes and the latter as the plasmon-like modes. Indeed, the 
latter ones (TE polarized) resemble the TM polarized surface plasmon polaritons of a 
thin metal film or a gap plasmon guide. We correlate the angular location of the peaks 
in resonant transmission to the corresponding modes in a bare NIM guide. In all our 
calculations, along with retention of both the sources of dispersion, we present results 
where material dispersion is suppressed. We show that such a step can lead to wrong 
conclusions as regards the nature and the magnitude of the delay. In fact, there is a 
danger of predicting superluminal transit while actually it is subluminal. It is believed 
that NIMs have to be essentially lossy [221 123]. In order to have some quantitative idea 
about the effect of absorption on the delay, we assume a causal response of the NIM 
medium. We show that even in presence of a much improved magnetic response it may 
be a difficult job to have high-Q guided and surface modes in the optical and near IR 
domains. Our study thus reiterates the important observation of Veselago in the context 
of absorption and dispersion of a realistic metamaterial. 

The structure of the paper is as follows. In section 2, we consider a generic negative 
index material slab in order to reveal the roles of intrinsic and structural dispersion. We 
show that there can be qualitative differences in the group delay if the material dispersion 
is suppressed. In section 3, we pick the NIM parameters from a recent experiment |21j, 
albeit with the assumption of small damping. We study group delay through a structure 
supporting resonant tunnelling. After a brief discussion of the NIM dispersion and the 
resulting delay in bulk, we investigate the plasmon-like modes of the structure. We show 
that the excitation of such modes can lead to enhanced delay as compared to FP modes, 
albeit at the expense of lower transmission. In section 4, we discuss the guided modes of 
the same structure, except for lower values of refractive index of the spacer layers. Like 
in the case of plasmon-like modes, these guided modes are also shown to lead to enhanced 
delay. In section 5, we fit the experimental data [21] for the magnetic resonance to a 
causal response and vary the relevant parameters in order to obtain better features of 
the NIM. We show that even a substantial improvement of the magnetic response fails to 
excite true high-Q guided and surface modes in the NIM guide. Finally, we summarize 
the main results in Conclusions. We make an important observation on how the large 
damping in the present day metamaterials can destroy the resonant tunnelling features, 
stressing the dire need to fabricate 'transparent' metamaterials. 

2. Group delay in a FP cavity containing a generic NIM 

Consider first the propagation of an optical pulse through a bulk NIM sample. Let the 
permittivity e and the permeability /i of the NIM be given by the following expressions 



HDH2H 




uj 2 — u)n + iVu 
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Recall that the group delay for a segment of length d in the bulk is given by 
d ( d\ , „ dn 



V 9 



n g , n g 



n{u) + uj—, (2) 



where n g and n are the group index and the refractive index of the material, respectively. 
It can be easily seen from (j2j) that group index n g is the same as the normalized delay 
rj (d/c). Since the source of delay in ([2j) is the material dispersion, we will refer to this as 
delay due to intrinsic dispersion. From a different viewpoint, delay r through a segment 
of length d can be linked to the frequency derivative of the phase <p t (accumulated over 
distance d) of the (amplitude) transmission coefficient t (= |t|e 4I?!,t ) [25] as 



T 



(3) 



where uo c is the carrier frequency of the pulse, r in ([3]) is also referred to as the Wigner 
phase time. 

As mentioned earlier, for finite structures, there is another source of frequency 
dependence of important parameters like transmission and reflection through it. For 
example, if one considers a slab of width d of a material with some refractive index, 
embedded in a medium with another index, the characteristic signatures of the Airy 
resonances of the FP cavity will be imprinted on the delay features. We will refer 
to such features (with neglect of material dispersion), as delay due to pure structural 
dispersion. The overall delay emerges as an interplay and competition of these two 
distinct sources of dispersion. It is thus clear that none of the sources of dispersion can 
be ignored for correct assessment of the delay. Irrespective of the nature of the delay, it 
can be evaluated by calculating the phase of the amplitude transmission coefficient and 
using the expression for the Wigner phase time. 

The group delay of a pulse passing through a NIM FP cavity was studied in 
detail and an analytical expression for the delay was derived [10]. It was shown that 
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Figure 1. (a) Real (solid) and imaginary (dotted) parts of the refractive index of the 
NIM, (b) the normalized delay for a bulk sample. The parameters chosen are: f p = 12 
GHz, f b = 6 GHz, f = 4 GHz and T/u = 1CT 3 . 
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Figure 2. (a) Intensity transmission coefficient T (top) and the normalized delay 
r/(d/c) (bottom) for a FP cavity as functions of frequency; (b) the same as functions 
of angle of incidence 9i for / = 5 GHz, T/llIq = 1CP 3 and d = 10 cm. The bottom 
right panel compares the normalized delays with (solid) and without (dashed) material 
dispersion. 

such a system can lead to large group delays at frequencies corresponding to the Airy 
resonances. A closer inspection of the expression of the Wigner phase time reveals 
that the delay has a very pronounced dependence on the frequency derivatives of the 
permittivity and the permeability. It is thus expected that negligence of the material 
dispersion would lead to quite erroneous conclusions, even to the extent of the sign of 
the delay. In other words, one may wrongly predict superluminal transit, while it is 
actually subluminal. 

In order to have a quantitative assessment of the contributions of intrinsic and 
structural dispersion for the delay, we first study the dispersive properties of the bulk 
material given by ([1]). Parameters were taken from the work of Smith et al |24| . namely 
f p = 12 GHz, f b = 6 GHz, /„ = 4 GHz and T/u = 1(T 3 , where f pJbfi = u Pibfl /{2ir). 
The frequency dependence of the real and imaginary parts of the refractive index is 
shown in figure UK In the bottom panel (figure [lb) we have plotted the time taken for 
propagation through a distance d in the bulk sample. It is clear that over the specified 
range of frequencies the bulk material exhibits normal dispersion which leads to a delay 
of the pulse. If propagation is through a slab of such a material, embedded in air, the FP 
resonances are imprinted on this background leading to large delays at these resonances. 
Now we demonstrate how suppression of the material dispersion can lead to completely 
erroneous conclusions. Consider the same FP slab, but now with frequency independent 
e = e/ =5 GHz, A* = / i /=5GHz- For reference we have chosen the material parameters at an 
intermediate frequency, namely, / = 5 GHz. The Airy resonances of this FP cavity 
are shown in the top panel of figure [2^. The bottom panel shows the corresponding 
normalized delay. As can be seen from this figure, just the structural dispersion of the 
slab for the said parameters leads to negative delay characteristics, which is just the 
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opposite of the actual situation. The same results can be verified from the angle scan 
of the transmission and the normalized delay at a particular frequency of incident light. 
In figure [2b we show the dependence of the intensity transmission T and the normalized 
phase time r/{d/c) as functions of the angle of incidence 0j at / = 5 GHz. For each 
angle of incidence we evaluated the complex amplitude transmission coefficients at two 
infinitesimally close frequencies to ± Auj and extracted the corresponding phases. The 
ratio of the differences in the phases to that of the frequencies lead to the phase time at 
us. Results for the cases with (solid line) or without (dashed line) intrinsic dispersion are 
shown in the lower panel of figure [2b . Material dispersion was suppressed by assuming 
e(oj ± Auj) = e(uj), fi(ou ± Auj) = fi(u>). Again it is clear that neglect of dispersion can 
lead to the incorrect result of negative delay. 



3. Delay in resonant tunnelling mediated by plasmon-like modes 

In this section we consider a structure which can lead to resonant tunnelling via the 
excitation of the surface (guided) modes. The structure is assumed to be symmetric 
(e» = 6f, ex — e 3 , d\ = d 3 ) with the central layer made up of NIM (see figure [3]). In 
order to establish a proper connection with the current experimental activities, and also 
to test their feasibility with our current proposal, we pick the NIM parameters like e(ui) 
and h{uj) from the experiment of Dolling et al [21] • As mentioned earlier, despite the 
huge volume of literature on novel metamaterials, very few give the complete dispersion 
data (real and imaginary parts of both e and /x) for the NIM [201 [23 EE]- Most of 
these data are extracted from transmission studies. The work of Dolling et al reported 
a figure of merit of about 3 at A ~ 1.4/xm, where Re(n) = —1, much needed for 
perfect lensing applications. We used the data digitized from their results for e and 
H with two significant changes. It is understood that a thin layer of (120nm in their 
experiment) metamaterial is highly anisotropic. In our calculations we assumed the NIM 
material to be isotropic and homogeneous. Hopefully 3d near isotropic metamaterials 
will be realized in near future. The second approximation concerns the losses in 
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Figure 3. Schematics of the layered structure with central NIM layer sandwiched 
between two spacer layers and high index prisms. All other materials except for the 
NIM are assumed to be non-magnetic. 
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the metamaterial. We later comment on how the actual losses of currently available 
metamaterials can wash out all the interesting effects reported here. Thus the other 
issue concerns the dire need of truly low-loss metamaterials. Instead of the actual 
experimentally observed losses, we assume a small loss of 27 at the working frequency 
ZJ for both permittivity and permeability. Thus at the specified frequency for the NIM 
we write 

e 2 (p) = Re(e d (uJ))+ij, 
/i 2 (uj) = Re(fi d (p)) + ij, 

where e^tU) and /id(uJ) are taken from the work of Dolling et al [2T] . The motivation 
for introducing the small damping is to understand, at least qualitatively, the effect 
of damping on the resonant tunnelling through a NIM guide. Henceforth, in this and 
the next section, we consider only the low-loss NIM guide with properties given by (j4j) 
excited by a TE-polarized light. 

In order to have the plasmon like modes of the NIM guide we choose the spacer 
layer refractive index n\ larger than \Re[ji2)\ at the working wavelength. In particular, 
we choose the following system parameters: A = 2ttuJ/c = 1.425/im,ej = 6.145, e\ = 
63 = 2.25, d\ = g?3 = lyum, d 2 = 2fim. The thickness of the spacer layers was chosen 
so as to optimize the excitation of a given mode of a NIM guide. The results for the 
intensity transmission (T), phase of transmission (<f) t ) and corresponding normalized 
delay (t/(c?t/c), dx = d\ + d 2 + d^) as functions of the angle of incidence are shown 
in figure @] a, b and c, respectively. The dashed, dotted and solid curves are for three 
different values of 7, namely, 7 = 0.0001, 0.001 and 0.01, respectively. The nature of 
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Figure 4. (a) Intensity transmission coefficient T (b) phase of transmission (in units 
of it) and (c) normalized delay as a function of the angle of incidence. The parameters 
used are A = 1.425/im,ei = e/ = 6.145, e\ = £3 = 2.25, di = ^3 = 1/^to, d-i = 2/im 
for 7 = 0.0001 (dashed), 0.001 (solid) and 0.01 (dotted). (Figure SJi inset) Plot of 
(/c| — k\)d 2 as a function of a^d as in the work of Wu et al [15] . (Figure 2J; inset) 
Normalized delays with (solid lines, values reduced by a factor of 700) and without 
(dashed lines) material dispersion for 7 = 0.001. 



Competition between structural and intrinsic dispersion in delay through a left handed mediumS 

the peaks in figure @] a and c can be easily understood if one recalls the critical angles 
for the relevant interfaces at the specified wavelength. For the prism/silica interface the 
critical angle is about 37.24° while that for the prism/NIM interface it is about 35.25°. 
For Oi > 37.24° waves are evanescent in silica as well as the NIM slab. Thus the peak 
occurring at 9 ~ 39° can be associated with the surface-plasmon-like mode of the NIM 
guide. Such a mode is the dual of the TM-polarized surface plasmon polariton in a 
metal film [19J and have been discussed in detail in [12], [19] and used even for QED 
applications [T5] . 

Following Wu et al [19] we now determine the symmetry of the mode occuring at 
9i = 39.45°. For this purpose we consider a symmetric NIM guide with width d 2 = 2fim 
embedded in silica (with e = 2.25). For such a guide with ^1/^2 — —0.7487 the 
symmetric (cosh-type) and anti-symmetric (sinh-type) mode dispersion curves are shown 
in the inset of figure UK- From the parameters used in figure EK, {k\ ~ kf)d 2 = —15.7594 
(horizontal line in the inset of figure UK). We obtain two values of a 2 d — (k 2 — k\) l l 2 d 
corresponding to the symmetric (solid curve) and anti-symmetric (dashed curve) modes. 
For the symmetric mode a 2 d ~ 5.949 (vertical line in the inset of figure UK). The 
excitation angle corresponding to this value of a 2 d is 39.65° and is closer to the value of 
9i from the RT data. Thus the resonance at 9i = 39.45° in figure UK is recognized as the 
symmetric plasmon-like mode. The peak close to 9i = 34.5° corresponds to one of the 
FP modes of the layered medium, since waves are propagating both in silica as well as 
in the NIM layer. It is clear from figure Hb that the phase of the transmission coefficient 
undergoes qualitatively different (positive and negative) changes as one sweeps the angles 
through the FP or the plasmon-like resonances. The oscillatory (plasmonic) modes are 
associated with positive (negative) jumps in phase. However, both lead to large delays. 
It is also clear how increasing damping gradually erases the RT features. For values of 
NIM losses as in [21], the transmission peak due to the 'plasmonic' mode gets washed 
completely though the enhanced delay properties survive. Of course, there is no use of a 
delayed signal if its peak amplitude is vanishingly small. FP modes and the associated 
delay still survive to some extent even for large losses. One other important aspect that 
should be noted from figure H] is that the plasmon-like modes can lead to much larger 
delays than the FP modes, due to the large quality factors associated with such modes. 

4. Delay in resonant tunnelling mediated by guided modes 

In this section, we consider again the RT configuration (see figure [3]), albeit with two 
changes. We choose a spacer layer with e x = e 3 = 1.0, so that n x < Re{n 2 ) and the 
guided modes could be excited. In other words, the spacer layer is chosen to be optically 
rarer than the NIM. The other change is in the width of the spacer layers, namely, 
d\ = d 3 = 0.25/im, in order to optimize the excitation of the modes. The critical 
angles for the prism/air interface is 23.79° and prism/NIM interface (at A = 1.425/zm) 
is 35.25°. For angles greater than 23.79° and less than 35.25°, the waves are evanescent 
in the spacer layers, while they are propagating in the NIM layer. Thus guided modes 
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Figure 5. Same as in figure 0] except that ei = 63 = 1, /Lti = fj,^ = 1, di = (I3 = 
0.25/im. Inset in figure[5^L shows the effective index of the modes as a function of guide 
width. Inset in figure [SJ; shows the normalized delays with (solid lines, values reduced 
by a factor of 200) and without (dashed lines) material dispersion for 7 = 0.001. 



can be excited in the NIM layer only for this range of angles. Any peak observed in 
the transmission profile in this range can be recognized as due to the excitation of these 
guided modes. 

We show the intensity transmission T, phase of transmission (f) t and the normalized 
delay rj (c?t/c) as functions of the angle of incidence 9i in figuresEh, b and c, respectively. 
The dashed, solid and dotted curves correspond to 7 = 0.0001, 0.001 and 0.01, 
respectively. It is clear from figure [5] that two guided modes corresponding to 0j = 
27.76° and 33.51°, respectively, are supported in the NIM layer. In the inset of figure 
[5^ we plot the effective index, n eS = k x /k , of the modes of a lossless NIM guide in 
air as a function of the guide width. Different curves are labeled by the corresponding 
mode numbers. Recall that for a NIM guide the m = mode does not exist and the 
lowest allowed mode (m = 1) has a cutoff (also noted in [T2J, [19] ). From the inset 
we confirm that for di = 2fim (vertical dashed line in the inset), only two guided 
modes (marked with circles) are supported. The excitation angles (9i) for these values 
of n e ff are 31.61° (for m = 1) and 27.28° (for m = 2), respectively, and they match 
well with the RT data. The mismatch can be attributed to the fact that losses were 
ignored. The loading of the guide with the prism in the RT configuration is another 
important source of this mismatch. The plot of the phase of transmission (figure [5b) 
shows positive jumps at the location of the guided mode resonances (compare with 
figure 0b). The plot of normalized delay (figure Eb) shows positive delays for both the 
guided modes. The intensity of transmission gets considerably weakened with increased 
damping though the delay characteristic is retained. In the inset of figure [5b we plot 
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the normalized delays for 7 = 0.001 with (solid line, values reduced by a factor of 200) 
and without (dashed line) including the effect of the material dispersion. We note again 
that consideration of material dispersion removes the erroneous negative delay features 
(obtained by disregarding the material dispersion) of NIM. Finally, a comparison of 
figures [5b and Hb reveal that the delay with the plasmonic mode is usually higher than 
that with the guided mode. This is perhaps due to the tighter confinement of the field 
and the larger quality factor of the plasmonic mode. 

5. Effects of losses in the framework of a causal response 

The treatment in the previous sections (3 and 4) at a given frequency was simple 
and losses were introduced just to demonstrate their devastating effects on resonant 
tunnelling features. The realistic treatment of losses over the whole spectral domain 
is much more complicated. In order to have a causal response, one must have 
susceptibilities analytic in the upper half complex plane, which is a direct consequence 
of the fact that response functions vanish for negative arguments. Moreover, the real 
and imaginary parts of the susceptibility are related by the Kramers-Kronig relation. 
Recent analysis of causality in the context of NIM dispersion sheds much insight on the 
underlying phenomena (with or without additional gain) J22J [23j [29] . These studies are 
based on the analytic properties of n 2 {uj) = e(u)(j,(u). A compensation or reduction 




Figure 6. The real (solid) and imaginary (dashed) parts of refractive index (top row), 
the normalized delay (middle row) and FOM (bottom row) as a function of A. The left 
top panel shows the fitted curves for digitized data from Dolling et al [5T]. The right 
panel shows the corresponding curves for improved magnetic resonance (/ = 10). 

of material losses narrows down the regime over which the left-handed behaviour is 
observed. On the other hand, as noted in the preceeding sections, increasing loss 
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washes out the interesting features associated with NIM guides. Attainment of a lossless 
NIM thus reaches an impasse. There should thus be a trade-off between the material 
loss and the regime of negative index behaviour for practical metamaterials. Magnetic 
losses are known to play a dominant role in controlling the negative index behaviour 
of metamaterials. In this section we analyze a material with Lorentz type magnetic 
permeability given by [30] 



fj,(u) = Hoo 5 ~T~T~- ' ( 5 ) 

where fj,^, u>o, F and 7 m have their usual meanings. Note that the experimental data 
of Dolling et al [21] can be fitted to (jSJ) with the following values of the parameters 
/ioo = 0.6, A = 2ttc/u =1A59 jum, 7 m /^ = 0.028 and F = 4.425 (j m /u) ). We used 
the digitized data for the permittivity. The dependence for the said parameter values 
is shown in the top panel of figure [6^,, while the middle and the bottom panels show 
the corresponding normalized delay in bulk sample and the FOM, respectively. It was 
mentioned earlier that in a guide with such material parameters, the RT features are 
completely washed out. Hence we probe the case of an NIM with improved magnetic 
response, whereby we scale down the damping 7 m by a factor f 1 with simultaneous 
scaling up of the oscillator strength F by f Q . The corresponding improved response for 
fo = f-v = 10 is shown in figure [6b. 
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Figure 7. Intensity transmission coefficient T (top row) , phase of transmission <j>t in 
units of 7r (middle row) and the time delays (bottom row) as functions of the angle of 
incidence 9i for the RT structure supporting plasmon-like modes (a) and guided modes 
(b). Parameters chosen for panel (a) are t\ = £3 = 4, d\ = d-j = 1/ito, g?2 = 0.5^tm and 
for panel (b) are t\ = £3 = 1, d\ = d% = 0.05/^m, c?2 = lfim. The common parameters 
are e, = e f = 6.145 and / =1 (solid), 10 (dashed), 100 (dotted). 
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We studied the transmission and pulse delay through a RT geometry for f Q — f y — 
10, 100 and compared it with the material of Dolling et al [21] (f — f 1 — 1). Delay in 
RT mediated by plasmon-like modes were studied in a structure with d% = d$ — Ifim, 
€\ = e 3 = 4 and d 2 = 0.5/im at A = 1.425/im and the results are shown in figure [Th. 
The critical angles for the prism/NIM interface is 44.42° and the prism/spacer layer 
interface is 53.79°. We show that plasmon-like modes can be supported in this structure 
with large time delays for improved magnetic responses. These features disappear for the 
experimental data of Dolling et al [2T] . A study of time delays mediated by guided modes 
with parameters di = d 3 = 0.05/im, ei = 63 = 1 and d 2 = 1.0/xm at A = 1.425/im (figure 
[7b) shows that the supported guided mode shows enhanced transmission and delay 
for improved magnetic responses (f D — / 7 = 10, 100) whereas the data corresponding 
to Dolling et al [21] (/ G = f 1 — 1) barely shows any transmission or delay features. 
Thus it should be possible to have enhanced delay characteristics with materials having 
improved magnetic responses. 

6. Conclusions 

We studied pulse delay through a layered medium containing a NIM layer. We 
considered three specific cases, namely, a generic NIM, a realistic, recently reported 
NIM with the assumption of low losses and NIM with improved magnetic response. We 
also assumed the NIM to be homogeneous. We showed that large delays in transmission, 
mediated by the modes of the NIM guide, are achievable exploiting a resonant tunneling 
geometry only for very low loss structures. Both guided and plasmon-like modes were 
exploited to this end. We probed the role of structural and intrinsic NIM dispersion 
in this delay and stressed the need to retain both for its correct assessment. We also 
commented on how the large damping of contemporary metamaterials can destroy the 
RT transmission features. We probed a NIM with better magnetic response to show the 
reemergence of these features. These calculations are carried out in the framework of a 
full causal theory. We thus ingeminated Veselago's observation about the importance 
and necessity of retaining intrinsic NIM dispersion in all relevant problems. We also 
emphasized the need for truly low-loss metamaterials for the realization of such compact 
delay devices. 

The authors are thankful to the Department of Science and Technology, 
Government of India, for supporting this work. SDG also thanks Girish S Agarwal 
for interesting discussions. 
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